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Abstract 

We present the first measurement of photoproduction of J ftp and of two-photon 
production of high-mass e+e~-pairs in electromagnetic (or ultra-peripheral) nucleus- 
nucleus interactions, using Au-|-Au data at ^/s^ = 200 GeV. The events are tagged 
with forward neutrons emitted following Coulomb excitation of one or both Au* 
nuclei. The event sample consists of 28 events with mg+g- > 2 GeV/c^ with zero 
like-sign background. The measured cross sections at midrapidity of da/dy {J/i^ + 
Xn,y = 0) = 76 ± 33 (stat) ± 11 (syst) /xb and dFu/dmdy (e+e" -I- Xn,y = 0) = 
86 ± 23 (stat) ± 16 (syst) ^b/(GcV/c2) for nZg+g- G [2.0,2.8] GeV/c^ are consistent 
with various theoretical predictions. 
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1 Introduction 



The idea to use the strong electromagnetic fields present in high-energy nucleus- 
nucleus collisions to study photoproduction at hadron colliders has attracted 
growing interest in recent years, see [1,2,3] for reviews. Electromagnetic inter- 
actions can be studied without background from hadronic processes in ultra- 
peripheral collisions (UPC) without nuclear overlap, i.e. impact parameters 
larger than the sum of the nuclear radii. This study focuses on the measure- 
ment of exclusively produced high-mass e'^e~ pairs at midrapidity in Au-|-Au 
collisions at ^/s^ = 200 GeV, Au + Au Au + Au + e"'"e~. The results have 
been obtained with the PHENIX detector [4] at the BNL Relativistic Heavy 
Ion Colhder (RHIC). 

The electromagnetic field of a relativistic particle can be represented by a 
spectrum of equivalent photons. This is the Weizsacker- Williams method of 
virtual quanta [5,6]. The number of photons in the spectrum is proportional 
to Z^, where Z is the charge of the particle, and the equivalent two-photon 
luminosity is thus proportional to Z^. The strong dependence on Z favours the 
use of heavy ions for studying two-photon and photonuclear processes. The 
virtualities of the equivalent photons when the field couples coherently to the 
entire nucleus are restricted by the nuclear form factor to = (u;^/7^ + qj_) < 
{Ii/RaY- Here, uo and g_L are the photon energy and transverse momentum, 
respectively, Ra is the nuclear radius and 7 the Lorentz factor of the beam. 
At RHIC energies, 7 = 108 and the maximum photon energy in the center- 
of-mass system is of the order of a;niax ~ 3 GeV corresponding to maximum 
photon-nucleon and two-photon center-of-mass energies of W^^^ ~ 34 GeV 
and - 6 GeV. 

The exclusive production of an e^e" pair can proceed either through a purely 
electromagnetic process (a two-photon interaction to leading order) or through 
coherent photonuclear production of a vector meson, which decays into an elec- 
tron pair. Exclusive photoproduction of vector mesons is usually thought of as 
proceeding via Pomeron-exchange, the perturbative-QCD equivalent of which 
is the exchange of two gluons or a gluon ladder. The Feynman diagrams for the 
two leading order processes are shown in Figure 1. The two-gluon picture is 
applicable to production of heavy vector mesons, such as the J/ip^ and to pro- 
duction of lighter mesons at high momentum transfers [7]. The J /ip production 
cross section is consequently a good probe of the proton [8] and nuclear gluon 
distribution, G^(a;, Q^), as well as of vector- meson dynamics in nuclear mat- 
ter [9,10]. For J/'^-production, the coverage of the PHENIX central tracking 
arm, —0.35 < rj < 0.35 corresponds to a range in the photon-nucleon center- 
of-mass energy between 21 < W^n < 30 GeV, with a mean (W^n) = 24 GeV. 
This corresponds to photon energies in the rest frame of the target nucleus of 
240 <E^< 480 GeV, with {E^) = 300 GeV. Mid-rapidity photoproduction of 
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Fig. 1. Lowest order Feynman diagrams for exclusive photoproduction of (a) J/ if; 
and (b) dielectrons, in ultra-peripheral Au+Au collisions. The photons to the right 
of the dashed line are soft photons that may excite the nuclei but do not lead to 
particle production in the central rapidity region. Both diagrams contain at least 
one photon and occur when the nuclei are separated by impact parameters larger 
than the sum of the nuclear radii. 

J/ip probes nuclear Bjorken-a; values of a; = fn^ji^/W^j^ ~ 1.5-10^^ [10], where 
the nuclear gluon density is partially depleted by "shadowing" effects [11] com- 
pared to the proton. 

The strong fields associated with heavy ions at high energies lead to large 
probabilities for exchanging additional soft photons in the same event. Most of 
these photons have too low energy to produce particles, but they can excite the 
interacting nuclei. The dominating excitation is to a Giant-Dipole Resonance 
(GDR) with energies 0(10 MeV), which decays by emitting neutrons at very 
forward rapidities, providing a very useful means to trigger on UPCs with Zero- 
Dcgrcc Calorimeters (ZDC). The probability for having a Coulomb excitation 
leading to emission of neutrons in at least one direction in coincidence with 
coherent J/iIj production is 55% ± 6% [12]. The probabilities for exchanging 
one or several photons factorise, i.e. the Coulomb tagging does not introduce 
any bias in the extraction of exclusive J/ ^l) photoproduction cross sections from 
these events [13]. The soft photons leading to moderate nuclear excitation are 
indicated to the right of the dashed line in Figure 1. Incoherent (or quasi- 
elastic) vector-meson photoproduction can also proceed via the interaction of 
the exchanged photon with a single nucleon in the nucleus. In that case, J/ if) 
photoproduction is always accompanied by nuclear breakup and emission of 
nucleons in the forward direction [14]. 

Photoproduction of vector mesons has been studied with lepton beams first in 
the 60s [15,16] and more recently at the electron- proton collider HERA [17,18]. 
Measurements of photonuclear production of p mesons [19,20], as well as of 
two-photon production of low-mass e+e~ pairs [21] in heavy ion interactions 
have been performed by the STAR collaboration. The CDF collaboration has 
studied two-photon production of e+e~ pairs [22] and exclusive production 
of n'^ ii~ pairs [23] in pp collisions at the Tevatron. The PHENIX analysis 
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presented here is the first on heavy final states in ultra-peripheral nucleus- 
nucleus collisions. It supersedes a preliminary study presented earlier [24]. The 
cross section for J/ip and e+e^ photoproduction are compared with various 
theoretical calculations [10,14,25,26,27]. 



2 Experimental setup 

The data presented here were collected with the PHENIX detector at RHIC 
during the 2004 high-luminosity Au+Au run at ^/s^ = 200 GeV. The PHENIX 
detector [4] , is a versatile detector designed to study the properties of strongly 
interacting matter at extreme temperatures and energy densities present in 
central heavy ion collisions. The current analysis demonstrates its capabili- 
ties to also study ultra-peripheral collisions, which have a very different event 
topology. The PHENIX central tracking system [28] consists of two arms, each 
covering \r]\ < 0.35 and A0 = 7r/2, equipped with multi-layer drift chambers 
(DC) followed by multi-wire proportional chambers (PC) with pixel-pad read- 
out. The tracking arms also have Ring-Imaging-Cerenkov (RICH, with CO2 
gas radiator) detectors [29] and electromagnetic calorimeters (EMCal) [30] for 
electron and positron identification. The PHENIX EMCal consists of six sec- 
tors of lead-scintillator sandwich calorimeter (PbSc, 15552 individual towers 
with 5.54 cm x 5.54 cm x 37.5 cm, I8X0) and two sectors of lead-glass 
Ccrenkov calorimeter (PbGl, 9216 modules with 4 cm x 4 cm x 40 cm, 
14.4X0), at a radial distance of ~5m from the beam line. 

The ultra-peripheral Au-|-Au events were tagged by neutron detection at small 
forward angles in the ZDC. The ZDC's [31,32] are hadronic calorimeters placed 
18 m up- and down-stream of the interaction point that measure the energy 
of the neutrons coming from the Au* Coulomb dissociation with ~20% energy 
resolution and cover |^| < 2 mrad, which is a very forward region^ . 

The events used in this analysis were collected with the UPC trigger set up 
for the first time in PHENIX during the 2004 run with the following charac- 
teristics: 

(1) A veto on coincident signals in both Beam-Beam Counters (BBC, cov- 
ering 3.0 < \ri\ < 3.9 and full azimuth) selects exclusive-type events 
characterised by a large rapidity gap on either side of the central arm. 

(2) The EMCal-Trigger (ERT) with a 2x2 tile threshold at 0.8 GeV. The 
trigger is set if the analog sum of the energy deposit in a 2x2 tile of 
calorimeter towers is above threshold (0.8 GeV). The efficiency for trig- 

^ Much larger than the crossing angle of Au beams at the PHENIX interaction 
point (0.2 mrad). 
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gcring at least one of the two high-energy e coming from the pair is 
estimated to be e'^^^ = 0.9 ± 0.1. 
(3) At least 30 GeV energy deposited in one or both of the ZDCs is required to 
select Au+Au events with forward neutron emission {Xn) from the (single 
or double) Au* decay. As has been discussed, this requirement leaves 
about 55% of the coherent and RilOO% of the incoherent J/ip events. 

The trigger efficiency is estimated to be String = 0.9 ± 0.1 from the ERT 
requirement. The inefficiencies due to requirements (1) and (3) were calculated 
and found to be negligible. 

The total number of events collected by the UPC trigger was 8.5 M, of which 
6.7 M satisfied standard data quality assurance criteria. The useable event 
sample corresponds to an integrated luminosity = 141±12 /xb~^ computed 
from the minimum bias triggered events. 



3 Data Analysis 

Charged particle tracking in the PHENIX central arms is based on a combi- 
natorial Hough transform in the track bend plane (perpendicular to the beam 
direction). The polar angle is determined from the position of the track in the 
PC outside the DC and the reconstructed position of the collision vertex [33]. 
For central collisions, the collision vertex is reconstructed from timing infor- 
mation from the BBC and/or ZDC. This does not work for UPC events, which, 
by definition, do not have BBC coincidences and often do not have ZDC co- 
incidences. The event vertex was instead reconstructed from the position of 
the PC hits and EMCal clusters associated with the tracks in the event. This 
gave an event vertex resolution in the longitudinal direction of 1 cm. Track 
momenta are measured with a resolution 6p/p ^ 0.7% © 1.0% p[GeV/c] in 
minimum bias Au-|-Au nuclear collisions [34]. Only a negligible reduction in 
the resolution is expected in this analysis because of the different vertex res- 
olution. 

The following global cuts were applied to enhance the sample of genuine 7- 
induced events: 

(1) A standard offline vertex cut \vtXz\ < 30 cm was required to select col- 
hsions well centered in the fiducial area of the central detectors and to 
avoid tracks close to the magnet poles. 

(2) Only events with two charged particles were analyzed. This is a restrictive 



For instance, for the low-pT photoproduced J/tp, the e decays have E « 
0.5 mj/^ 1.6 GeV. 
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criterion imposed to cleanly select "exclusive" processes characterised 
by only two isolated particles (electrons) in the final state. It allows to 
suppress the contamination of non-UPC (mainly beam-gas and peripheral 
nuclear) reactions that fired the UPC trigger, whereas the signal loss is 
small (less than 5%). 

Unlike the J/ip ^ e+e" analyses in nuclear Au+Au reactions [34,35] which 
have to deal with large particle multiplicities, we did not need to apply very 
strict electron identification cuts in the clean UPC environment. Instead, the 
following RICH- and EMCal-based offline cuts were used: 

(1) RICH multiplicity >2 selects which fire 2 or more tubes around 
the track within the nominal ring radius. 

(2) Candidate tracks with an associated EMCal cluster with dead or noisy 
towers within a 2x2 tile are excluded. 

(3) At least one of the tracks in the pair is required to pass an EMCal cluster 
energy cut {Ei > 1 GeV \\ E2 > 1 GeV) to select candidate in the 
plateau region above the turn-on curve of the ERT trigger (which has a 
0.8 GeV threshold). 

Beyond those global or single-track cuts, an additional "coherent" identifi- 
cation cut was applied by selecting only those e~^e~ candidates detected in 
opposite arms. Such a cut aims at reducing the high-px pairs while improving 
the detection of the low-px pairs expected for 77, 77! production. Neverthe- 
less, after all the previous cuts were applied the infiuence of this selection is 
found to be small; there is only one event in which the e"^ and e~ are in the 
same arm and have mg+g- > 2 GeV/c^. 

Finally, J/ip were reconstructed by invariant mass analysis of the measured 
e'^e" pairs. There was no remaining like-sign background after the aforemen- 
tioned analysis cuts. 

The cross sections were obtained after correcting the raw number of signal 
counts for the geometrical acceptance of our detector system, and the efficiency 
losses introduced by the previously described analysis cuts. Acceptance and 
efficiency corrections were obtained using the PHENIX GEANtS [36] simula- 
tion package with input distributions from the STARLIGHT Monte Carlo (MC), 
based on the models presented in [12,25,37]. The measured j+p V+p cross 
sections from HERA and fixed target experiments with lepton beams are used 
as input to the models, starlight well reproduces the existing N / dyd(l)dpT 
distribution of coherent p production in UPC Au-|-Au events measured at 
RHIC by STAR [19,20]. Helicity conservation is assumed in the model, and 
the angular distribution of the decay products {J/ip e+e^) is given by 
dN/dcos{9) oc 1 -|- cos^(6') in the J/ip center-of-mass. The angular distribu- 
tion is different from that for p production followed by the decay p — > tt'^tt". 
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because of the different spin of the daughters, as well as from the angular 
distribution in two-photon interactions 7 + 7 e^e~ . We generated 5-10"^ co- 
herent J/ if) and 8-10^ coherent high-mass e^e~ pairs (mg+e- > 1 GeV/c^) in 
Au-|-Au collisions accompanied by forward neutron emission. The simulated 
events were passed through the same reconstruction programme as the real 
data. 

Table 1 summarises the J/ip and dielectron acceptance and efficiency cor- 
rection factors obtained from our simulation studies. For instance, for J/-^ 
photoproduction the correction is 1/(2.49 ±0.25)%, of which the experimental 
acceptance to detect the decay electron pair is about 5% (for J/ip produced 
at \y\ < 0.35). In the 77 e^e~ sample, most of the electrons/positrons are 
emitted at very forward angles. The fraction of events with Impair | < 0.35 and 
2.0 < rUg+e- < 2.8 GeV/c^, where both the electron and positron are within 
\r]\ < 0.35 is 1.10%. The corresponding numbers for 2.0 < nie+e- < 2.3 GeV/c^ 
and 2.3 < m^+e- < 2.8 GeV/c^ are 1.11% and 1.08%, respectively. The accep- 
tance and efficiency corrections have a systematic uncertainty of 10% resulting 
from the accuracy of the simulation to describe the detector, the electron iden- 
tification parameters, and the event vertex position resolution. 

Table 1 

Coherent J/i/j and e+e~ (continuum) acceptance and efficiency for |ypair| < 0.35 as 
a function of invariant mass raiig(\ Tlic^ las! lii u^ sliows IIk^ ( rigger efficiency. 



me+e- [GeV/c2] 


Acc X e 




J/^P 


(2.49 ± 0.25) • 


10-2 


e+e- [2.0,2.8] 


(2.24 ±0.22) • 


10-3 


e+e- [2.0,2.3] 


(2.16 ±0.22) • 


10-3 


e+e- [2.3,2.8] 


(2.33 ± 0.23) • 


10-3 


(■trigg 


0.9 ±0.1 





4 Results and Discussion 

The measured e+e^ invariant mass distribution for the sample is shown in 
Figure 2 a). The amount of background can be estimated from the number 
of like-sign events (i.e. events where two electrons or two positrons are recon- 
structed). We find no like-sign pairs for me±g± > 2 GeV/c^, compared with 
28 events with an e^e~ pair with rrie+e- > 2 GeV/c^. The shape is consis- 
tent with the expected contribution from the two processes in Figure 1: a 
continuum distribution corresponding to two-photon production of e+e" pairs 
and a contribution from J/t/j ^ e+e". Since the offline cuts {Ei > 1 GeV || 
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a) b) 




iTiee (GeV/c^) me3(GeV/c=) 



Fig. 2. Left: a) Invariant mass distribution of e~^e~ pairs fitted to the combination 
of a dielectron continuum (exponential distribution) and a J ftp (Gaussian) signal. 
The two additional dashed curves indicate the fit results with the maximum and 
minimum continuum contributions considered in this analysis (see text), b) J/ip 
invariant mass distribution after subtracting the fitted dielectron continuum signal 
in a). 

E2 > 1 GeV) cause a sharp drop in the efficiency for rrie+e- < 2 GeV/c^, we 
include only pairs with m^+^- > 2 GeV/c^ in the analysis. 

The invariant mass distribution is fitted with a continuum (exponential) curve 
combined with a Gaussian function at the J/ip peak, as shown by the solid 
curve in Figure 2 a). Simulations based on events generated by the starlight 
MC (see last paragraphs of Section 3) processed through geant have shown 
that the shape of the measured continuum contribution is well described by an 
exponential function dN/dm^+e- = A ■ e'^"^<=+<=' . Those simulations allow us to 
fix the exponential slope parameter to c = —1.9 ±0.1 GeV~^c^. The combined 
data fit is done with three free parameters: the exponential normalisation {A), 
the J/ ip yield and the J/ ip peak width (the Gaussian peak position has been 
fixed at the known J /ip mass of mj/^ = 3.097 GeV/c^ [38]). Figure 2 b) shows 
the resulting invariant mass distribution obtained by subtracting the fitted 
exponential curve of the dielectron continuum from the total experimental 
e^e~ pairs distribution. There is a clear J/ip peak, the width of which (cxj/^ ~ 
155 MeV/c^) is consistent with the J/ip width from our full MC. 

The J /ip and continuum yields and the corresponding statistical errors are cal- 
culated from the fit. The dashed curves in Figure 2 a) show the maximum and 
minimum e^e~ continuum contributions considered, including both the statis- 
tical and systematical uncertainties. The systematic uncertainties were deter- 
mined varying the dielectron continuum subtraction method using a power-law 
form instead of an exponential function and by modifying the corresponding 
fitted slope parameters by ±3(T. The propagated uncertainty of the extracted 
yields was estimated to be one count in both cases. The total number of J/ip^s 
is: Nj/^ = 9.9 ± 4.1 (stat) ± 1.0 (syst), and the number of e+e" continuum 
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Table 2 

J /ip ^ and e+e" continuum yields obtained from the fit of the data to an 

exponential plus Gaussian function per invariant mass range. Systematic errors are 
obtained as described in the text. 

rrie+e- [GeV/c^] Yield 







= 9.9 ±4.1 (stat) ± 1.0 (syst) 


e+e- [2.0,2.8] 




= 13.7 ± 3.7 (stat) ± 1.0 (syst) 


e+e" [2.0,2.3] 




= 7.4 ± 2.7 (stat) ± 1.0 (syst) 


e+e- [2.3,2.8] 




= 6.2 ± 2.5 (stat) ± 1.0 (syst) 



pairs for m^+^- G [2.0,2.8] GeV/c^ is: 7Ve+e- = 13.7 ± 3.7 (stat) ± 1.0 (syst). 
Table 2 shows the obtained results per invariant mass range. 

Figure 3 a) shows a scatter plot of invariant mass mg+g- vs. pair px- From the 
plot, it is clear that most of the pairs outside the J j'^ peak originate in coher- 
ent processes with very low pair transverse momenta (pt ^ 100 MeV/c), as 
expected for two-photon interactions. For events with mg+e- around the Jjil) 
mass, however, there are a few counts at larger values which cannot be as- 
cribed to background events, since there are no like-sign pairs above 2 GeV/c^. 
A purely coherent production - corresponding to events where the fields cou- 
ple coherently to all nucleons and the nucleus remains in its ground state 
(7 -I- A — > V ^ A) - would yield % 200 MeV/c after reconstruction. On 
the other hand, incoherent production (7 + A V ^ X) - dominated by 
the quasi-elastic vector meson production off one nucleon inside the nucleus, 
7-1-iV— >-y + iV - results in much larger px for the photoproduced J j'^ [14]. 
The cross sections for coherent and incoherent J/V" photoproduction in UPCs 
at RHIC are expected to be of the same order [14] . We discuss below whether 
our data confirm such a prediction. 

The transverse momentum distribution of the events with mg+e- > 2 GeV/c^ 
is shown in Figure 3 b). For clarity, only points below px < 1 GeV/c are drawn. 
The px is here the magnitude of the vector sum of the px of the electron and 
positron. One sees a clear enhancement of events with very low transverse 
momenta, consistent with coherent production. The squared form-factor of a 
gold nucleus. 



sin(i?px) — i?px cos(-Rpx) 
(i?Px)Ml + (apx)2) 



is shown for comparison. Here, R = 6.7 fm is the gold radius, and a = 0.7 fm 
represents the difTuseness of the nuclear surface [39]. The magnitude of the 
form-factor is a free parameter fitted to reproduce the spectra. Figure 3 also 
presents the corresponding distribution expressed in terms of the squared mo- 
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A = 1282.4 ± 256.8 (x';NDF= 41.3/99) 
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M e [2.7,3.5] GeV/c= 
dHJdp', = A*|F(p^)|= 
A = 20.4± 5.9 (x^/NDF= 17.5/99) 
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p2 (GeV^/c^) 



Fig. 3. Top: (a) Scatter plot of e+e" mg+g- vs. pair px- (b) dN/dpT distribution 
of the pairs with mg+g- > 2 GeV/c^ fitted to the Au nuclear form factor, Eq. (1). 
Bottom: dN/dp^ distribution of pairs with (c) rrig+g- G [2.0,2.7] GeV/c^ and (d) 
rrig+g- G [2.7,3.5] GeV/c^, fitted to the expected Au nuclear form factor. Note the 
difference in scale on the x-axis in the lower two plots. 

mentum transfer from the target nucleus, ~ ~'^-> events with mg+g- 
corresponding to the dielectron continuum mg+g- G [2.0,2.8] GeV/c^ (Figure 3 
c) and below the J/'?/'-peak, mg+g- G [2.7,3.5] GeV/c^ (Figure 3 d). 

The extracted yields of J/ip and e"'"e~ are used to calculate the final cross sec- 
tion for photoproduction at midrapidity in ultra-peripheral Au+Au collisions 
accompanied by forward neutron emission. For dielectrons at midrapidity {y 
is the rapidity of the pair) the double differential cross section is: 



e+e-+Xn 



dy drrie 



1 



1 



[s/|<0.35,Am^ 



Acc ■ e ■ Etriqa ' C.int Ay Amg+g 



86 ± 23 (stat) ± 16 (syst) /ib/(GeV/c^) for mg+g- G [2.0,2.8] GeV/c^ 



For J/'?/' at midrapidity the differential cross section is: 



dy 



|j/|<0.35 



1 

'br 



N 



j\.CC ' E ' Sifiqq ' H 



•'int 



1 



76 ± 31 (stat) ± 15 (syst) /xb. 



(3) 
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The correction factors (and corresponding uncertainties) are quoted in Table 1 
as described in previous sections, and BR = 5.94% is the known J/ip dielec- 
tron branching ratio [38] . Table 3 summarises the measured cross sections per 
invariant mass interval. 

Table 3 

Measured J/il^ and e^e" continuum photoproduction cross sections at midrapidity 
in ultra-peripheral Au+Au collisions (accompanied with forward neutron emission) 
at ^/s^ = 200 GeV. The rightmost column in the lower part shows the starlight 
predictions [37]. 



da/dy\y=o [/jh] 




76 ± 31 (stat) ± 15 (syst) 


me+e- [GeV/c2] 


(fa/dmf,+f,-dy\y=o [//b/(GeV/c^)] 




data STARLIGHT 


e'^e~ continuum [2.0,2.8] 


86 ± 23 (stat) ± 16 (syst) 90 


e+e~ continuum [2.0,2.3] 


129 ± 47 (stat) ± 28 (syst) 138 


e+e- continuum [2.3.2.(S] 


60 ± 21 (slat) ± 11 (syst) 61 



The measured dielectron cross sections at midrapidity are in very good agree- 
ment with the STARLIGHT predictions for coherent dielectron photoproduc- 
tion (rightmost column of Table 3) [37]. Exclusive dilepton production in 
STARLIGHT is Calculated combining the two equivalent (Weizsacker-WiUiams) 
photon fluxes from each ion with the Breit- Wheeler formula for 77 — > 
The agreement between STARLIGHT and other leading order calculations [40] 
is good as long as the pair invariant mass is not too low. A recent calcula- 
tion has found that higher order terms suppress the e'^e~ cross section by 
29% in the invariant mass range 140 < mg+e- < 165 MeV/c^ [41]. A re- 
duction of the same magnitude in the invariant mass range considered here, 
2.0 < rUe+e- < 2.8 GeV/c^, would still be in agreement with our measurement. 

The final J/ip + Xn cross section is compared to the theoretical predictions 
computed in references [12,14,25,27,37,42] in Figure 4. The rapidity distri- 
butions of Strikman et al. [14] and Kopeliovich et al. [27] have been scaled 
down according to [12] to account for the reduction of the yield expected 
when requiring coincident forward neutron emission {Xn). The scaling has 
been applied as a function of rapidity with the integrated cross section being 
55% of the original one. The band covered by the Strikman et al. predictions 
includes the J/'0 cross sections with and without gluon shadowing (as imple- 
mented in the Frankfurt-Guzey-Strikman, FGS, prescription [43]). Strikman 
and Kopeliovich predictions for the coherent and incoherent photoproduc- 
tion cross sections are drawn separately in Figure 4 a) and summed up in 
Figure 4 b). STARLIGHT [37] and Gongalves-Machado [42] calculations only 
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evaluate the coherent contribution. 



As mentioned above, the measured pair pt distributions suggest coherent J/V' 
photoproduction (7 + ^4 — J/ip + X) and a possible additional incoherent 
(7 + A?" — > J/'^-\-X) contribution at higher px- To give an indicative estimate 
of the size of the incoherent contribution, we can assume that it corresponds 
to the counts in the J/il) mass window with p\ > 0.1 (0.05) GeV^/c^. This 
corresponds to about 4 (6) counts, which amounts to a contribution of about 
40 (60)% of the total J/i/j production, compatible with the theoretical calcu- 
lations [14]. The limited data statistics prevents us from separating in a more 
quantitative way the two components. Note that although the acceptance cor- 
rection for the J/ip was calculated using a Monte Carlo which includes only the 
coherent component, the obtained correction is also a reasonable approxima- 
tion for the incoherent component, provided that quasi-elastic scattering on a 
single nucleon, j + N ^ V + N, gives the main contribution. The polarisation 
of the vector meson will then be the same as for coherent production, and the 
reduction in acceptance because of the different px range will be of the order 
of ~10-20%. If the incoherent contribution to the total J/ip photoproduction 
was 40%, the coherent J/ip cross section would become ~ 46 /ih. 

Despite these uncertainties, the final J/ip cross section is in good agreement, 
within the (still large) statistical errors, with the theoretical values computed 
in [12,14,25,27,37,42] as shown in Figure 4. The current uncertainties unfortu- 
nately preclude any more detailed conclusion at this point regarding the two 
crucial ingredients of the models (nuclear gluon shadowing and J/ip nuclear 
absorption cross section). The statistical uncertainties can be improved with 
significantly higher Au+Au luminosities and a concurrent measurement of the 
J/ ip in the dimuon decay channel in the more forward acceptances covered by 
the PHENIX muon spectrometers, as collected in year 2007 and expected in 
the future. 

Finally, one can attempt to compare the obtained photonuclear J/ip cross 
sections to those from e-p collisions at HERA by dividing the measured dif- 
ferential cross section {da/dy) with the (theoretical) equivalent photon spec- 
trum (dN^/duj). At midrapidity: a^^^j^^^ = {da^^^j^^^^/dy)/{2dN^/dw), 
with 2N^ — 6.7 (10.5) for the coherent (incoherent) spectrum at a photon- 
nucleon center-of-mass energy of (VF^p) = 24 GeV. Assuming, for the sake 
of simplicity, a 50% - 50% contribution from coherent and incoherent inter- 
actions in our total ultra-peripheral J/ip sample, the extracted photonuclear 
cross sections are: a{'y+Au J/t/j+Au) — 5.7±2.3 (stat) ±1.2 (syst) /xb, and 
(7(7 + Au^ J/i> + X) ^ 3.6 ±1.4 (stat) ± 0.7 (syst) /ib, respectively. A fit to 
the results from the HI and ZEUS collaborations [17,18] over their measured 
energy range gives (7(7 + p ^ J/ip +p) = 30.5 ±2.7 nb at W.yp = 24 GeV. 
Therefore, the ratios a{-f + Au J/^)/a{'^ + p J/^) = 186 ±88, 118 ±54 
for the coherent and incoherent components (statistical and systematic errors 
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assumed independent and added in quadrature) are consistent with a scal- 
ing of the photonuclear cross section with the number of nucleons in gold 
{A = 179): a(7 + Au ^ J/^/j) = A'' ai^y + p ^ J/ip) with a^oh = 1-01 ±0.07, 
and aincoh = 0.92 ±0.08, respectively^ . 



UPC Au-Au ^ J/'F Xn 200 GeV 

• PHENIX 

Starlight, coh. 

Strikman et al, coh. 

Strikman et ai, incoh. 

— ' — ' Goncalves-iVlachado, coh. 

— Kopeiiovich et ai, coh. 

■ iii-iii-iii-ii Kopeiiovich et ai, incoh. 



1:1: 

II 
It 



I 



-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

a) y 



T5 



UPC Au-Au ^ J/1' Xn 200 GeV 



PHENIX 

Starlight, coh. 

Strikman et ai, coh. + incoh. 

Goncalves-iVtachado, coh. 

Kopeiiovich et ai, coh. + incoh. 




Fig. 4. Cross section of J/tp + Xn production at midrapidity in UPC 
Au+Au collisions at ^s^^ = 200 GeV compared to theoretical calculations: 
STARLIGHT [12,25,37], Strikman, et al. [14], Goncalves-Machado [42], and Kope- 
iiovich, et al., [27]. The error bar (box) shows the statistical (systematical) uncer- 
tainties of the measurement. When available, the theoretical calculations for the 
coherent and incoherent components are shown separately in a), and summed up in 
b). 



5 Summary and Conclusions 



We have presented the first exclusive photoproduction of J/ip e^e^ and 
high-mass e^e~ pairs ever measured in nucleus-nucleus (as well as hadron- 
hadron) interactions. The measurement has been carried out by the PHENIX 
experiment in ultra-peripheral Au±Au interactions at ^fs^ = 200 GeV tagged 
by forward (ZDC) neutron detection from the (single or double) Au* dissocia- 
tion. Clear signals of J/ ip and high mass dielectron continuum have been found 
in the data. We have observed 28 e~^e~ pairs in rrie+e- £ [2.0,6.0] GeV/c^ with 
zero like-sign background. Their spectrum is peaked at low ~ 90 MeV/c 

^ Note, for comparison, that repeating the same exercise for the photoproduced p in 
the STAR UPC measurement [20], a{'y+Au p+Au) = 530±19 (stat)±57 (syst) ^b 
for (W-yN) ~ 12.5 GeV, and taking the experimentally-derived value ct(7 + p — > 
p + p) = 9.88 ph from [25], yields Ucoh = 0.75 ± 0.02 closer to the A^/^-scaling 
expected for soft particle production. 
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as expected for coherent photoproduction with a reahstic Au nuclear form fac- 
tor. 



The measured number of continuum e'^e events in the PHENIX acceptance 
for me+e- e [2.0,2.8] GeV/c^ is: N^e+e') = 13.7 ± 3.7 (stat) ± 1.0 (syst). After 
correcting for acceptance and efficiency losses and normalising by the measured 
luminosity, we obtain a cross section of (Pa /drrie+e-dy {e^ e~ -\- Xn)\y=Q = 86 ± 
23 (stat) ± 16 (syst) iih/ (GeV/c^), which is in good agreement with theoretical 
expectations for coherent exclusive dielectron production in photon-photon 
interactions. 



The measured invariant mass distribution has a clear peak at the J/iIj mass 
with an experimental width in good agreement with a full GEANT-based sim- 
ulation for UPC production and reconstruction in the PHENIX detector. The 
measured number of J/ilj mesons in the PHENIX acceptance is: N(J/il)) — 
9.9 ±4.1 (stat) ±1.0 (syst). The higher distribution suggests an addi- 
tional incoherent contribution to J/ip photoproduction in accordance with 
predictions [14], but statistical limitations prevent a more quantitative esti- 
mate. After correcting for acceptance and efficiency losses and normalising 
by the measured luminosity, the total J/ip photoproduction cross section is 
da/dy{J/ip + Xn)\y=o = 76 ± 31 (stat) ± 15 (syst) fih, which is consistent 
(within uncertainties) with theoretical expectations. The low background in 
the present data sample shows that future higher luminosity runs with reduced 
experimental uncertainties of the measured cross sections will provide more 
quantitative information on the nuclear gluon distribution and J/ip absorption 
in cold nuclear matter at RHIC energies. 
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